Uhlmann fidelity between two-mode Gaussian states 
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We analyze the Uhlmann fidehty of a pair of n-mode Gaussian states of the quantum radiation 
field. This quantity is shown to be the product of an exponential function depending on the relative 
average displacement and a factor fully determined by the symplectic spectrum of the covariance 
matrix of a specific Gaussian state. However, it is difficult to handle our general formula unless the 
Gaussian states commute or at least one of them is pure. On the contrary, in the simplest cases 
n = 1 and n = 2, it leads to explicit analytic formulae. Our main result is a calculable expression 
of the fidelity of two arbitrary two-mode Gaussian states. This can be applied to build reliable 
measures of quantum correlations between modes in various branches of quantum physics. 
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I. INTRODUCTION. 

One of the most important issues in quantum informa- 
tion is the closeness between input and output states of 
various quantum protocols When the input state is 
pure a natural figure of merit of any protocol is the input- 
output quantum-mechanical transition probability. The 
situation is more complicated in the mixed-state case. 
Fortunately, an appropriate description in such a case 
could be connected with a concept introduced long ago 
by Uhlmann 0. Let us consider two arbitrary states, p' 
and /5", of a given quantum system. When both states 
are mixed, Uhlmann asserted that a suitable measure 
of the likeness between their properties is the maximal 
quantum-mechanical transition probability between their 
purifications in an enlarged Hilbert space 2, 3]. This 
is an extended notion of transition probability between 
quantum states which is now called fidelity 3] and has 
the intrinsic expression Q 



F{p',p") 



Tr 



(1.1) 



Indeed, when at least one of the quantum states is pure, 
the fidelity (|l.ip reduces to the familiar overlap between 
the states: T{p',p") = Tr(/5'p''). Important properties 
listed and discussed in Refs.[2|-|8| highly recommended 
fidelity as a measure of distinguishability of quantum 
states. Especially useful is the property put forward in 
Ref. [3| that fidelity equals the minimal squared overlap 
of the probability distributions for the outcomes of any 
general measurement. Fidelity-based metrics such as the 
Bures metric 0, Q and related ones [1-0| have proven 
to be fruitful in quantum optics and quantum informa- 
tion. For instance, the Bures metric provides insightful 
distance-type degrees of entanglement [lol [ll| , nonclas- 
sicality [12], and polarization [13|t3- Moreover, in the 
framework of the theory of parameter estimation, Braun- 
stein and Caves [13] established a distinguishability met- 
ric which is proportional to the Bures distance between 



neighboring quantum states. Their findings in Ref.[l7j 
successfully apply in recent developments of quantum es- 
timation theory, as shown in the review [l8l | . Besides be- 
ing an unavoidable ingredient in quantum information, 
fidelity is exploited in various areas of physics where the 
concept of quantum transition probability plays a key 
role. For instance, in condensed matter physics, fidelity 
is frequently employed as an indicator of criticality in 
quantum phase transitions [l^, [2^ . 

In spite of playing a key role in quantum mechanics 
due to its operational meaning, fidelity between mixed 
states has a rather limited use so far. This happens be- 
cause, in general, Uhlmann's expression (jl.ip is not easy 
to calculate even on finite-dimensional Hilbert spaces 
Within continuous- variable settings, explicit calcu- 
lations of fidelity were performed for one-mode Gaussian 
states (GSs) in Refs.[2l|, 111], for some special two-mode 
GSs in Refs.[ll|, or for commuting two- mode states in 
Refs.0,[l^. 

In what follows, we concentrate on the GSs of 
continuous-variable systems, which are important for 
both quantum optics and quantum information process- 
ing [23]. From the experimental side, they are quite ac- 
cessible and were effectively obtained with light, Bose- 
Einstein condensates, trapped ions, and Josephson junc- 
tions. On the other hand, the GSs of the quantum ra- 
diation field are important resources in several quantuirr 
information protocols [13, HI]. In particular, from the 
theoretical point of view, the two- mode GSs are a perfect 
test bed for studying entanglement [Tli o r other kinds of 
correlations between the two modes [2y, [S^] . Therefore, 
an exact analytic formula for the fidelity of an arbitrary 
pair of two-mode GSs is highly desirable and its deriva- 
tion constitutes a main goal of the present paper. Let us 
mention that we have previously derived the fidelity be- 
tween two-mode symmetric and squeezed thermal states 
and used it in evaluating a Gaussian degree of eirtangle- 
ment defined with the Bures metric JJj]. Nevertheless, 
an explicit formula of the fidelity between arbitrary two- 
mode GSs was not yet found. We acknowledge, however. 
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Ref. [28| as a useful prospect for the multimode Gaussian 
case. 

In order to study the fidehty of two GSs, we find it 
convenient to introduce the operators: 



A 



'■I -11 
PP , 



B 



WW, 



C:=v^V^, (1.2) 



with the properties B = C'^C and Tv{A) = Tr{B) ^ 0. 
When this trace is strictly positive, Tr(^) > 0, then a 

state 



PB ■■= [Tr(6)]-ie, 
exists, so that the fidelity (II.ip reads: 

J-(p',p") = Tr(i) [Tr(V^) 



(1.3) 



(1.4) 



Equation \1A\ is the starting point of our explicit evalu- 
ation of the fidelity between two-mode GSs. 

From now on, we assume that p' and p" are multimode 
GSs. The plan of the paper is as follows. In Sec. II we 
recall the phase-space description of the GSs and men- 
tion Williamson's theorem, which is at the heart of our 
method. Some useful properties of the product of two 
Hilbert-Schmidt Gaussian operators are collected in Sec. 
III. Based on the composition rules of such a product, 
we then proceed in Sec. IV to derive the fidelity (jl.4|) of 
a pair of multimode GSs. We express the second factor 
in the formula (11.41) in terms of the symplectic eigenval- 
ues of the covariance matrix (CM) of the GS pB- We 
develop our treatment in the multimode framework as 
far as possible. Then we restrict ourselves to the work- 
able case of two commuting multimode GSs. In Sec. V, 
the general expression of the fidelity of two GSs enables 
us to retrieve a previous explicit formula for one-mode 
states, in agreement with Refs.[23, Then we reach 
the main goal of this work: we provide a versatile and 
workable formula for the fidelity between two-mode GSs. 
An interesting application of the two-mode formula con- 
cludes Sec. V. The final section deals with a discussion 
of our analytic results and their potential use in studying 
quantum correlations between modes. 



II. PHASE-SPACE DESCRIPTION OF THE 
GAUSSIAN STATES. 

An n-mode Hilbert-Schmidt operator F is defined by 
its Weyl expansion: 



1 



(27r)'^ 



d^"w f{u) D{-u), (2.1) 



where u is a vector in the Euclidian space M^" , the weight 
function f{u) is square-integrable over R^", and D{u) is 
an n-mode Weyl displacement operator: 



D{u) :-(g)i^, ((«),). 



(2.2) 



Let qj and Pj be the canonical quadrature operators of 
the jth field mode, and qj and pj their eigenvalues. In 
Eq. (j2.2|) . {u)j denotes a vector in R^, whose components 
are precisely qj and pj : 



iu)j := 



U2J-1 
U2j 



(j = l,2,...,n). (2.3) 



Further, Dj{{u)j) is the corresponding single-mode Weyl 
displacement operator, 

Dj{{u)j) := exp [~i {qjp, - pjq.j)] . (2.4) 

Recall the standard matrix J of the symplectic form 



on 



p2n 



which is block-diagonal and skew-symmetric: 



J := Jfc, Jk := 



k=l 



1 

-1 



(A; = l,2,, 



(2.5) 



Two properties of the Weyl operators (|2.2p are useful: 
• composition law: 



D{u) D{v) = exp --u^ Jv D{u + v) 



orthonormalization relation: 



(2.6) 



Tr 



b\u)b{v) \ = (2^)"5(2n)(^ _ (2.7) 



For instance, the latter yields the scalar product of two 
Hilbert-Schmidt n-mode operators ()2.1|) . 



Tr 



1 



(27r)" 



(2., 



(2.9) 



as well as the trace-formula for a weight function, 



/(u) = Tr b{u)F 



Note the value /(O) = Tr(F) and that an operator ((2?T|) 
is self-adjoint iff f*{-u) = f{u). 

We focus on multimode Hilbert-Schmidt Gaussian op- 
erators (HSGOs). By definition, the weight function i 
of such an operator is a shifted Gaussian, 



1 



(2.10) 



with a complex shift vector ^ e C^" and a symmetric 
matrix J- G M2„(C) whose real part is positive definite: 
J- + T* > 0. This implies that the square matrix J- is 
invertible, so that the Gaussian weight function (I2.10p is 
integrable 



^2n 



u exp 



—iP'Tu — if'^u 
2 ^ 



-(2^)"[det(^)]-^exp(--e^ T-'C], 



3?e 



{[det(J-)]^}>0' 



(2.11) 
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A special example of HSGO is the density operator p 
of an n-mode GS, which is positive and of unit trace: 
Tr(/5) = 1. Its weight function (|2.9|) . called the charac- 
teristic function (CF) of the state, is fully determined by 
the first- and second-order moments of all the quadrature 
operators: 



X(u) = exp 



-]^u^Vu-i{{u)^Yu 



(2.12) 



In Eq. (|2TT2)) . V G M2„(R) denotes the symmetric CM 
of the GS p. It fulfills the concise Robertson-Schrodinger 
uncertainty relation, 



Ct(v + ^j)C^o, 



(CeC 



2n\ 



(2.13) 



as a necessary and sufficient condition [30|,|31[. Therefore, 
the CM V is positive definite, V > 0, and, according to 
Williamson's theorem [32| , is congruent via a symplectic 
matrix S G Sp{2n, M) to a diagonal matrix: 




det(V) 



e 



det 1 V-f - J 



n 



(2.14) 



In Eqs. (j2.14p . the positive numbers kj, {j — 1, 2, . . . , n), 
are the symplectic eigenvalues of the CM V [33|, and all 
the single-mode matrices Ij are equal to the 2x2 iden- 
tity matrix. As a consequence, the n-mode Robertson- 
Schrodinger uncertainty relation (j2.13p is equivalent to 
the inequalities Kj ^ i, (j = 1, 2, . . . , n). 



III. PRODUCT OF TWO HILBERT-SCHMIDT 
GAUSSIAN OPERATORS. 

We consider the product F = F'F" of two n-mode HS- 
GOs whose weight functions f'{u) and f"{u) are Gaus- 
sians, Eq. (PTH]) . Making use of Eqs. (g^, and 
(|2.11[) . one finds that the weight function (|2.9I) of the 
product F has the same form. Although we do not list 
the corresponding composition rules |34i] here, we just 
stress that the real part of the resulting symmetric ma- 
trix T G M2„(C) is positive definite in two special cases 
we are interested in: 

1. Both factors F' and F" are self-adjoint. 

2. The factors are adjoints to one another: F' = 
(F")t. 

Hence, in both above-mentioned cases, the product F is 
itself an n-mode HSGO. In the first case, F is not self- 
adjoint, unless the two factors commute. In the second 



one, is a positive operator, proportional to the density 
operator of a GS whose CM is precisely F G M2„(R). 

We further investigate the product A :~ p' p" of two 
n-mode GSs. First, making use of Eqs.(2.8)-(2.11), we 
get its trace: 



Tr(i) = [det (V' + V")]" 



X exp 



^-(<5(^.))^(V' + V")-'<5(^i) 



> 0. 



(3.1) 



We have introduced above the relative average n-mode 
displacement := (u)p' — (u)p". 

Other properties of the product operator A are found 
by evaluating its weight function ()2.9|) via Eq. (|2.ip used 
for both its factors p' and p" . After finally performing 
a routine integral of the type (|2.11l) we are left with 
the Gaussian weight function /a(u), Eq. ()2.10p . which is 
determined by the following composition rules: 

y^A = -\j + (v" + ^ j) (V + vy' (v + ^ j) ; 

U = (u)^, - (v' ~ ^ (V + Vr' S{u). (3.2) 

Similarly, we employ the Hilbert-Schmidt scalar product 
(|2.8p in conjunction with the Gaussian integral (|2.1ip in 
order to get the formula 



Tr(A^) = Tr (A) 2"" [det (FaT'' 



(3.3) 



An essential non-trivial step is that the determinant of 
the matrix Fa, Eq. p.2|) . always factors as follows: 



det (Fa) 



det [(JVQ {JV") - 1/] 
det (V + V") 



(3.4) 



The determinant (|3.4p easily emerges when arranging 
the matrix Fa in a product form. On the other hand, 
Eq. (|3.2p implies the identity 



J) det (V" I J) 



/ i \ dct(V' 

^^^H-^'+S-'r det(V' + V") 

(3.5) 

Both the above determinants (13. 4p and (|3.5I) are mani- 
festly symmetric with respect to the GSs p' and p" . 

We exploit the composition rules (|3.1|) and p.2|) in 
three special situations. First, when the GSs coincide, 
p' = p" =: then Fa = Vb Va- Hence, Eqs. ((XT|) 
and (13.21) read: 



Tr(p2) = 2-" [det{V)r^ ^ 1: 



Va = ^ ( V - i JV-i J 
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(3.6) 



(3.7) 



Second, Eqs. ([3J)) . ([3J| and ((2?T4)) enable us to write 
down two conditions, both of them necessary and suffi- 
cient in order that an n-mode GS a is pure. They are 
expressed in terms of its CM V, as follows: 



4 



Scalar condition: 



Tr ((7^) = 1 ^ dct(V) = 2-2" 
1 
2 



'«j = 77, (i = l,2,...,n); 



• Matrix condition: 

Va = V < 



(jvy 



I. 



(3.8) 



(3.9) 



Third, the square root \/p of a given n-mode GS p with 
the CF (I2.12P is proportional to an equally displaced n- 
mode GS: its CM V determines the proportionality factor 



Tr(V/5) = \/det(2V), 
and satisfies Eq. p.7p 

1 1 



(3.10) 



V = - V- - JV^V 



Accordingly, its symplectic eigenvalues determined 
by those of the CM V, hereafter denoted Kji 



^ ^^j + \l i^jf - \, (i = l,2,...,n). (3.12) 



IV. 



FIDELITY OF TWO MULTIMODE 
GAUSSIAN STATES. 



We successively apply Eqs. ((XTOl) . ((2lH) . and (I3TT2)) 
to the GS ps, Eq. (|1.3p . for getting a general expression 
of the fidelity (IlII) : 



:F{p',p")=Tt{p'p") 



[(«b) 



(4.1) 



In Eq. the set {(k_b)j , (j = 1, 2, . . . , n)} is the 

symplectic spectrum of the CM Vb of the n-mode GS 
Pb- Some expected properties of the fidelity can be read 
on the above expression. For instance, let us note the 
inequalities 



T{p',p")^Tv{p'p")>0, 



(4.2) 



with saturation of the first one when at least one of the 
GSs p' and p" is pure. The second inequality is displayed 
by Eq. the Hilbert-Schmidt scalar product of any 

pair of n-modc GSs is strictly positive. Another property 
of the fidelity (14. ip is its symplectic invariance: 

J- {u{S)p'uHS), U{S}p"U\S)) = J-(p', p"), (4.3) 

where S £ Sp{2n, M) and U (S) are the unitary operators 
of the metaplectic representation on the Hilbert space of 
the n-mode states. 



There is a special case when the multimode fidelity 
(|4.ip can be readily evaluated, namely, that of commut- 
ing GSs. As a prototype of such states, we mention the 
n-mode thermal states. We concentrate on a pair of com- 
muting n-mode GSs, p' and p", whose CMs V and V" 
have the symplectic spectra {k^} and {k"}, respectively. 

The commutation relation [p' , p"] = implies, via the 
resulting equality Vb = J^a, that the diagonalizable ma- 
trices JV" and JV" commute, so that they have a com- 
plete system of common eigenvectors. Since the matrix 
JVb has the same eigenvectors, we find the composition 
law of the symplectic eigenvalues: 



(j = l,2,. 



, ,nj 



(4.4) 



Accordingly, the fidelity formula (|4.ip simplifies to: 



(3.11) ^(/3',p")=n 




(4.5) 



Equation (|4.5p displays two additional properties of the 
fidelity, namely, its symmetry, T{p" ,p') = T[p\p"\ and 
the saturable inequality 



J-(p',p")^l: J-(p',p")-l iff 



(4.6) 



We come back to the general case of two non- 
commuting n-mode GSs: writing an explicit formula of 
their fidelity (|4.ip seems to be a hard task. Indeed, eval- 
uation of the CM Vb of the n-mode GS ps, Eq. 
requires a three-time application of the composition rule 
(1221) and therefore the CM Vb & M2„(M) has a compli- 
cated structure. However, we succeeded to calculate its 
symplectic invariants det(VB) and det (Vb + ^J)- In- 
deed, taking advantage once again of Eqs. p.Sp and 
(|2.1ip . we find a formula which is similar to Eq. 



Tr(62^ = Tr{B) ^ 2"" [det (Vs)]"^ 



(4.7) 



Their comparison gives, via Eq. (|1.2p . the expression 

On the other hand, a rather involved calculation makes 
use of Eqs. p.2p and (|3.1ip to yield the determinant 

^ \ det(V- + f j)det(V" + |j) 
det[VB + -Jj^ + V") ' 



(4.9) 



which coincides with the no n- negative determinant (13.51) . 

When at least one of the GSs p' and p", say p' , is pure, 
then insertion of the appropriate conditions (13.81) and 
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(13:9)) into Eq. gS]) leads to the equality det (Vb) = 2"^", 
showing that the GS ps is also pure. 

For later convenience, let us introduce the following 
notations: 



det (V + V") ^ 1; 



2^" det 



(JV') {JV")--I 



> A: 



A:=22"det( V' + ^Jjdctfv" + ^j) ^0. (4.10) 



The inequality A ^ 1 originates in the feature of the pos- 
itive definite matrix V := 5 (V + V") of being the CM 
of an n-mode GS determined up to its average displace- 
ment. Let us remark that our main results in this section 
are the explicit expressions of the symplectic invariants 
(|4.8p and (14. 9|) for GSs with any n. However, only for 
one- and two-mode states these invariants are sufficient 
to get analytic formulae for the fidelity via Eq. (|4.ip . 



V. EXPLICIT RESULTS 

A. One-mode Gaussian states. 

For the sake of completitude we begin with the one- 
mode case. The symplectic spectrum of the CM Vb con- 
sists of a single eigenvalue, 



According to Eqs. (|5^ and (|CTI)) . the fidelity (0|) 
reads: 



^{P:P") = exp 



n/aTa-Va 



(5.1) 



The symmetric formula (j5.ip was derived and analyzed 
gradually long ago in Refs. [22j. In addition, by using 
the structure of the one-mode GSs, which are displaced 
squeezed thermal states [11], we succeeded to prove ex- 
plicitly in the present work the property ()4.6|) of their 
fidelity (jSll) . 



B. Two-mode Gaussian states. 



mula 

F{p,p") 



exp 



\{5{u)f{V' + Vr'5{u) 



(5.2) 



Previous results reported in Refs.jll] for special states 
(squeezed thermal and symmetric) are seen to be partic- 
ular cases of Eq. (|5.2p . 

The determinants (|4.10p of 4 x 4 matrices occurring 
in Eq. (j5.2p are easily computable. As a salient example, 
we evaluate them for a pair of undisplaced two-mode GSs 
whose CMs V and V" are both in an unsealed standard 
form. This means that they are partitioned into 2x2 
diagonal submatrices of the type 



hi C 

C bol 



I = 



1 
1 



c = 



c 
d 



(5.3) 



Here hi ^ i, 62 = 5, c ^ A matrix (|5.3p has two 
symplectic invariants: 

det(V) = (6162 - 



det I V + -J 



c^) (6162 - d^) > 0, 

det(V)-i(6? + 62 + 2cd) + 1^0. 

(5.4) 

We denote b'i,b2,c',d' and b'{,b2,c",d" the standard- 
form parameters of the CMs V and V", respectively. 
Then, Eqs. (|ilI71) . (lO)) . and ((5^ allow us to factor 
both determinants A and F and to visualize the sym- 
plectic invariant A: 



A = 



{b[ + br)ib'2 + b'^)-ic' + c"y 



(b'l+b'l) {b'2 + b'^)-{d' + d'r 



T = 16■ 



{[b[b'2-{d'f] [b'lb'^-icJ'f 

lib[b'l + b'2b'^ + 2d'c") + l: 
4 16 



ib[b'( + b'^b'^ + 2c' d" 



1 

16 



A = 16 |det(V') - \ [{b'lf + [b'^f + 2c'd'] + ^1 
X |det(V") - i [(6'/)' + + 2cV'] + 1 



(5.5) 



We are left to substitute the formulae (|5.5p into Eq. ([5 
written for undisplaced two-mode GSs. 



The CM Vb has two symplectic eigenvalues, {hb)i and 
{1^3)2- Their squares are the roots of a monic quadratic 
binomial whose coefficients are expressed in terms of the 
determinants (|4.10l) . A straightforward calculation ex- 
ploiting Eqs. (|3.2p and (|4.ip leads to the symmetric for- 



VI. CONCLUSIONS 

To sum up, in this work we have tackled a long- 
standing problem in quantum information with contin- 
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uous variables: the Uhlmann fidelity of two multimode 
GSs. Our simple phase-space approach led to the im- 
portant formulae (|^- (|331) and (|4Jl) - (|49|) which are 
valid for any number of modes and could be useful for fu- 
ture research. However, we stress that an analytic result 
for 71-mode GSs seems hard to be found when 71 ^ 3 : 
indeed, one then needs to solve a higher-degree charac- 
teristic equation that cannot be easily obtained. 

While in the single-mode case an explicit formula was 
found many years ago, we derived here a computable an- 
alytic expression in the two- mode case, Eq. (|5.2p . We 
expect that this formula will open a productive research 
in quantifying bipartite quantum correlations (entangle- 
ment and discord) in the Gaussian scenario. An exam- 
ple at hand is the already performed evaluation of a de- 



gree of entanglement for some important classes of GSs 
(squeezed thermal states, symmetric states) in Refs.[ll|. 
The analysis of an arbitrary two-mode GS requires our 
present result. Use of fidelity in evaluating the Gaussian 
geometric discord defined with the Bures metric is cur- 
rently one of our main interests. Note also that fidelity 
in mixed-state estimation or reconstruction could be a 
potential important tool. It appears that applications of 
the fidelity to the emerging field of quantum information 
and to other interesting areas of quantum physics are far 
from being exhausted. 
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